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Abstract. We present recent progress in the calculation of the helium fine-structure splitting 
of the 2 3 Pj states, based on the quantum electrodynamic theory. Apart from the complete 
evaluation of ma 7 and m 2 /M a 6 corrections, we have performed extensive tests by comparison 
with all experimental results for light helium-like ions and with the known large nuclear charge 
asymptotics of individual corrections. Our theoretical predictions are still limited by the 
unknown ma s term, which is conservatively estimated to be 1.7 kHz. However, comparison 
with the latest experimental result for the 2 3 Po — 2 3 P2 transition [M. Smiciklas and T. Shiner, 
Phys. Rev. Lett. 105, 123001 (2010)] suggests that the higher-order contribution is in fact much 
smaller than the theoretical estimate. This means that the spectroscopic determination of a can 
be significantly improved if another measurement of the 2 3 Po — 2 3 P2 transition in helium-like 
Li+ or Be 2+ ion is performed. 



1. Introduction 

The quantum electrodynamic (QED) theory of atomic energy levels has achieved a precision 
level that makes possible the determination of nuclear properties, like the charge radius, the 
magnetic dipole, or even the nuclear polar izability from measured atomic spectra. If the nuclear 
structure effects are negligible or can be eliminated, one may obtain fundamental constants 
from comparison of theoretical predictions with experimental results. The most important 
examples include the Rydberg constant determined from hydrogen spectroscopy, the electron 
mass derived from the bound-electron g factor in hydrogen-like ions, and a obtained from the 
helium fine structure. As first pointed out by Schwartz in 1964 [I], the splitting of the 2 3 Pj 
levels in helium can be used for an accurate determination of the fine structure constant a. The 
attractive features of the fine structure in helium as compared to other atomic transitions are, 
first, the long lifetime of the metastable 2 3 Pj levels (roughly two orders of magnitude longer than 
that of the 2p state in hydrogen) and, second, the relative simplicity of the theory. Schwartz's 
suggestion stimulated a sequence of calculations [2-5] , which resulted in a theoretical description 
of the helium fine structure complete up to order ma® (or a 4 Ry) and a value of a accurate to 
0.9 ppm [6]. 

The present experimental precision for the fine-structure intervals in helium is sufficient for 
a determination of a with an accuracy of 14 ppb from Refs. 0[8] and even 5 ppb from Ref. [9]. 



In order to match this level of accuracy in the theoretical description of the fine structure, the 
complete calculation of the next-order, ma 7 contribution and an estimation of the higher-order 
effects is needed. Work towards this end started in the 1990s and extended over two decades 
[10-19]. In 2006 the first complete evaluation of the ma 7 correction to the helium fine structure 
was reported by one of us (KP) [20J. However, the numerical results presented there were in 
disagreement with the experimental values by more than 10 standard deviations (a). 

In our recent investigations [21|,I22|. we recalculated, using formulae from Ref. [20J, all effects 
up to order ma 7 to the fine structure of helium and performed calculations for helium-like ions 
with nuclear charges Z up to 10. The calculations were extensively checked by studying the 
hydrogenic {Z — >• oo) limit of individual corrections and by comparing them with the results 
known from the hydrogen theory. We found several problems in previous numerical calculations 
and, in the meantime, the experimental value of the 2 3 Pi — 2 3 P2 transition was changed by 3a 
[8]. As a result, the present theoretical predictions are in agreement with the latest experimental 
data for the fine-structure intervals in helium, as well as with most of the experimental data 
available for light helium-like ions. Our calculation of the ma 7 correction for the fine-structure 
splitting in light helium-like atoms was reported in Refs. \22\ I23j. In this paper, we present a 
detailed description of all corrections to helium fine structure and a summary of the numerical 
results. 



2. QED theory of the helium fine structure 

According to the quantum electrodynamic theory (QED) the energy levels of an atomic system 
are a function of the fine structure constant a and the electron-nucleus mass ratio. We omit 
possible nuclear structure effects, as their contribution to the helium fine structure is negligible. 
The fine-structure splitting Ef s (a) can be expanded in powers of a, 



n(5) 
v fs 



p(6) 



+ 4 7) + 0(a S 



(1) 



The expansion terms E, 



(n) 



m 



a n g(ri) are Q £ orc j er m a n . They implicitly depend on the 



electron-nucleus mass ratio and may additionally involve powers of In a. The advantage of this 
approach is that each of the expansion terms is expressed as the expectation value of some 
effective Hamiltonian, as presented in the following. For convenience, we first consider the 
infinite nuclear mass limit, and then account for the finite nuclear mass corrections separately. 

The dominant contribution to the helium fine structure is induced by the spin-dependent 
part of the Breit-Pauli Hamiltonian, which is, for an infinitely heavy nucleus, 
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where f = r\—r-i- The above Hamiltonian includes the effect of the anomalous magnetic moment 
(amm) a e , which is given by [23] (neglecting small vacuum-polarization corrections coming from 
particles heavier than an electron) 
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Expanding the amm prefactors in Eq. ([2]), Hf s can be written as a sum of operators contributing 
to different orders in a 
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whereas H^ amm and -£ff s 7 a mm are the amm parts of the corresponding higher-order operators. 



Here, H^ s and H^ s are the complete effective Hamiltonians to order m a 4 and m a 5 , respectively, 

whereas H^ mm and #£Lm are the amm 
The contributions to the fine structure are 

= (Hjp) + 0(m/M), (5) 
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where the expectation values are calculated with the corresponding eigenstate of the 
nonrelativistic Hamiltonian Hq 
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The finite nuclear mass corrections up to order ma 5 are conveniently divided into three parts, 
termed the mass scaling, the mass polarization, and the recoil operators. The effect of the mass 
scaling is accounted for by including the prefactor (m r /m) 3 in the operator Hf s , where m r is 
the reduced mass for the electron-nucleus system. The effect of the mass polarization can be 
accounted for to all orders by evaluating expectation values of all operators on the eigenfunctions 
of the Schrodinger Hamiltonian with the mass-polarization operator (m r /M)pi ■ p% included. 
The third effect is induced by the recoil addition to the Breit-Pauli Hamiltonian 
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3. The spin-dependent ma 6 contribution 

The m a 6 contribution to the helium fine structure is a sum of the second-order perturbation 
corrections induced by the Breit-Pauli Hamiltonian and the expectation value of the effective 
fine-structure Hamiltonian to this order, , 
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Here, 1/(E$ — Hq)' is the reduced Green function and HjQ is the spin-independent part of the 
Breit-Pauli Hamiltonian, 

n$ = -(pf+rfJ + ^^W + ^-^^ + ^ri, (io) 

where we have omitted a term with <5 3 (r) since it vanishes for the triplet states, ill consists 
of 15 operators first derived by Douglas and Kroll (DK) [2] in the framework of the Salpeter 
equation. These operators were later rederived using the much simpler effective field method in 
Rcf. [15 . The result is 
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where the Bi are given in Table 1. 

The finite nuclear mass corrections to the ma 6 contribution can be divided into the mass 
scaling, the mass polarization, and the operator parts. The mass scaling prefactor is (m r /M) 4 
for the B 2 , B3, B4, and B$, (m r /M) 5 for the other Bi operators, (m r /M) 6 for the second- 
order corrections involving the first term in Eq. (jlOp . and (m r /M) 5 for all other second-order 
corrections. The mass polarization effect is most easily accounted for by including the mass 
polarization operator in the zeroth-order Hamiltonian. The operator part comes from recoil 
corrections to H^\ H^, and H^\ The recoil part of is given by Eq. (JSj) - The spin- 
independent recoil part of the Breit-Pauli Hamiltonian is 

^-iiEK(f+f)«+ rf >. <*> 

Recoil corrections to the DK operators were studied by Zhang [T2| and by Pachucki and 



Sapirstein [19]. The result is given by the effective Hamiltonian 
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4. The spin-dependent ma 7 correction 

The ma 7 correction to the helium fine structure can be conveniently separated into four parts 
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The first term above combines all terms with InZ and In a [11-13,15,20], 
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The second part of £^ is induced by effective Hamiltonians to order ma 7 . They were derived 
by one of us (K.P.) in Refs. |20U21j . (The previous derivation of this correction by Zhang I12j 
turned out to be not entirely consistent.) The result is 
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The Hamiltonian Hq is induced by the two-photon exchange between the electrons, the electron 
self-energy and the vacuum polarization. It is given by [20] 
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Here, the terms with In Z compensate the logarithmic dependence implicitly present in the 
expectation values of the singular operators 1/r 3 and 1/r 5 , so that matrix elements of Hq do 
not have any logarithms in their 1/Z expansion. The singular operators are defined through 
their integrals with the arbitrary smooth function / 
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where 7 is the Euler constant. The effective Hamiltonian Hh represents the anomalous magnetic 
moment (amm) correction to the Douglas-Kroll ma 6 operators and is given by [20] 
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where the Hi are presented in Table [TJ The last term of £ i in Eq. (|16p , the Hamiltonian 
-^i 7 L m is the ma 7 amm correction to the Breit-Pauli Hamiltonian, see Eq. IHb. 
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The third part of £^ is given by the second-order matrix elements of the form [20] 
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where H^ is the effective Hamiltonian responsible for the nonlogarithmic ma correction to 
the energy 
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-f^fs is the spin-independent part of the Breit-Pauli Hamiltonian given by Eq. ()10p . and H^ is 
the ma 5 amm correction to H^\ see Eq. ([2|). 

The fourth part of £ ^ is the contribution induced by the emission and reabsorption of virtual 

(7) 

photons of low energy. It is denoted as £ L and interpreted as the relativistic correction to the 

(7) 

Bethe logarithm. The expression for £\ reads [16] 
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where 5 (. . .) denotes the first-order perturbation of the matrix element (. . .) by H^ \ implying 
perturbations of the reference-state wave function, the reference-state energy, and the electron 
Hamiltonian. 



5. Results for helium fine-structure 

Summary of the individual contributions to the fine-structure intervals of helium is given in 
Table Numerical results are presented for the large foi and the small vyi intervals, defined by 

vol = [E(2 3 P )-E(2 3 P 1 )]/h (24) 
u 12 = [E(2 3 P 1 )-E(2 3 P 2 )]/h. (25) 

We note that the style of breaking the total result into separate entries used in Table [2] differs 
from that used in the summary tables of the previous papers by Pachucki et al. \2Q\ I21j . In 
particular, the lower-order terms listed in Table III of Ref. [20] and in Table II of Ref. [21] 
contained contributions of higher orders, whereas in the present work the entries in Table [2] 
contain only the contributions of the order specified. 

A term-by-term comparison with the independent calculation by Drake was performed in 
Ref. [23J. We observe good agreement between the two calculations for the lower-order terms, 



Table 2. Summary of individual contributions to the fine-structure intervals in helium, in 
kHz. The parameters [25] are a' 1 = 137.035 999 679(94), cRoo = 3 289 841960 361(22) kHz, and 
m/M = 1.370 933 555 70 x 10~ 4 . The label (+m/M) indicates that the corresponding entry 
comprises both the non-recoil and recoil contributions of the specified order in a. 
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° Ref. [7J, b Ref. [26], c Ref. [22], d Ref. [8], e Ref. [28], f Ref. [9]. 

namely, for the ma , ma 5 , and ma 6 corrections. However, for the recoil correction to order 
ma 6 , our results differ from those of Drake by about 0.5 kHz for both intervals. The reason 
for this disagreement seems to be different for the large and the small intervals. For the large 
interval, the deviation is due to the recoil operator part, whereas for the small interval, it is 
mainly due to the mass polarization part (see discussion in Ref. [21]) . 

Our present estimates of the uncalculated higher-order effects for helium are larger than 
those in the previous studies [17} 118], The previous estimates were significantly less than 1 kHz. 
They were based on logarithmic contributions to order ma 8 corresponding to the hydrogen fine 
structure. However, a larger contribution might originate from the nonlogarithmic relativistic 
corrections. So our present estimate is obtained by multiplying the ma e contribution for the 
^02 = ^oi + ^12 interval by the factor of (Za) 2 , which yields a conservative estimate of ±1.7 kHz 
for all foil ^12) and ^02 intervals. All nuclear structure effects are completely negligible at the 
current precision level. The finite nuclear size correction is estimated to yield 18 Hz for foi and 
6 Hz for v\2- 

Our result for the vq\ interval of helium agrees well with all recent experimental values 
[3 [261 [27]. F° r the v\2 interval, theoretical result is by about 2a larger than the values obtained 
in Refs. J7[ EE] but in agreement with the latest measurement by Hessels and coworkers [8]. 
Our theoretical prediction for the vq2 interval is in excellent agreement with the very recent 
measurement of Smiciklas and Shiner [9j. Comparison with this experimental result suggests 
that the higher-order contribution might in fact be much smaller than our conservative estimate. 
This means that, if an independent measurement on Li + or Be 2+ confirms the smallness of the 
ma 8 terms, the helium determination of a will be significantly improved. The measurement 
should be performed for the 2 3 Pq — 2 3 P2 transition, since it is not affected by the singlet-triplet 
mixing effects, which strongly depend on Z. 

In summary, the theory of the fine structure of helium and light helium-like ions is now 



complete up to orders ma 7 and a 6 m 2 /M. The theoretical predictions agree with the latest 
experimental results for helium, as well as with most of the experimental data for light helium- 
like ions. A combination of the theoretical and experimental results [9] for the 2 3 Pq — 2 3 i-2 
interval in helium yields an independent determination of the fine structure constant a 

a' 1 = 137.035 999 55(64)(4)(368) , (26) 

where the first error is the experimental uncertainty, the second one is the numerical uncertainty, 
and the third comes from the estimate of the ma 8 term (±1.7 kHz). The result (|26p is accurate 
to 27 ppb and in agreement with the recent value obtained from the electron g factor [29) . 
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